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ABSTRACT 

An experimental apparatus to measure the particle density 
distribution in the freeboard of an atmospheric fluidized bed uas 
designed and constructed. The density versus height measured by 
the sampling apparatus gives a similar exponential decrease as 
previous investigations have found. 

A particle trajectory model is developed uhich calculates the 
height and particle density distributions above the bed surface of 
an atmospheric fluidized bed. The parameters input to the model 
are the superficial velocity, initial partical velocity, gas jet 
velocity and duration, and the particle size distribution of the 
bed mass. The model uas evaluated using the experimental data for 
jet velocity, duration, and particle size. The predicted slope of 
the particle density versus height in the freeboard agrees uith 
the experimentally measured slope uithin Z0%. 

A sensitvity analysis using the trajectory model, resulted in 
a determination of the particle distributions in the freeboard of 
a fluidized bed as affected by varying the input parameters to the 
trajectory model. The most significant effects uere achieved when 
the jet velocity or duration uas altered. 



Thesis Supervisor 
Title: 



Dr. Leon Glicksman 
Senior Research Scientist 



ACKNOWLEDGEMENTS 



This thesis was prepared under the guidence and supervision 
of Dr. Leon Glicksman. His encouragement , patience, advice and 
understanding uere all greatly appreciated, for without then, this 
uork would not have been possible. 

The help and support of Ton Yule, who freely gave his tine, 
expert advice and assistance, is nost strongly appreciated. I 
sincerely thank him for his efforts. 

A special word of appreciation is given to the Tennessee 
Valley Authority for funding this research and making it all 
possible. 

Most of all, I would like to thank my fiancee Heide, for her 
support, understanding, assistance and love during the entire 
project. She was always there when I needed her, even when her 
own work was laborsome. 



3 



A DETERMINATION OF PARTICLE DENSITY DISTRIBUTIONS 
ABOVE FLUIDIZED BeOS 

TABLE OF CONTENTS 

ABSTRACT „ . . Z 

ACKNOWLEDGEMENTS 3 

TABLE OF CONTENTS 4 

LIST OF FIGURES G 

LIST OF TABLES 12 

CHAPTER I - INTRODUCTION 14 

CHAPTER II - PARTICLE SAMPLING APPARATUS Z2 

Design Alternatives Z2 

Apparatus Requirements 24 

Apparatus Design 2S 

Apparatus Testing 37 

CHAPTER III - EXPERIMENTAL PROCEDURE 40 

Fluidized Bed Conf igurat i on 40 

Equipment Set-up 43 

Sampling Procedure 49 

Sample Analysis 53 

CHAPTER IV - COMPUTER MODEL SB 

Introduction SB 

Model Theory 56 

Testing of Program S3 

CHAPTER V - EXPERIMENTAL RESULTS AND CONCLUSIONS 68 

Minimum Fluidization Velocity 68 

Entrainment Analysis 68 

Particle Size Distribution 80 

Oscilloscope Trace Analysis 85 

Sample Ueignt vs Bed Activity Correlation 89* 



4 



CHAPTER VI - TRAJECTORY MODEL RESULTS AND DISCUSSION . 9G 



Selection of Baseline Parameters SB 

Typical Output Using Baseline Parameters 39 

Model Sensitivity Analysis 10G 

variation of superficial velocity Uo 108 

variation of initial particle velocity Upo ... 115 

variation of jet velocity Uj 121 

variation of jet duration tj 128 

variation of particle distribution in bed mass 134 
Comparison of Model uith Experimental Results ... 140 

CHAPTER VII - CONCLUSIONS AND RECOMMENDATIONS 145 

Conclusions 145 

Recommendations 14G 

REFERENCES 148 

APPENDIX A - MOMENT OF INERTIA CALCULATIONS FOR 

PADDLES 150 

APPENDIX B - ERROR DETERMINATION OF VACUMN COLLECTION 

SYSTEM 15G 

APPENDIX C - SAMPLE TRAP CLOSURE TIME TEST 158 

APPENDIX D - SOLENOID TORQUE AND DYNAMIC ANALYSIS ... 1G2 

APPENDIX E - PARTICLE SIZE DISTRIBUTION ANALYSIS 1S9 

APPENDIX F - MEAN BED FLOW VELOCITY DETERMINATION ... 171 

APPENDIX G - COMPUTER PROGRAM FOR MEAN BED FLOW 

VELOCITY CALCULATION 175 

APPENDIX H - PARTS LIST FOR APPARATUS 178 

APPENDIX I - LISTING OF ALL SAMPLE DATA OBTAINED 182 

APPENDIX J - IMAGE ANALYZER OUTPUT Z0Z 

APPENDIX K - OSCILLOSCOPE TRACES 218 

APPENDIX L - COMPUTER PROGRAM FOR TRAJECTORY MODEL .. 235 

APPENDIX M - ANEMOMETER CALIBRATION 248 



5 



LIST OF FIGURES 



Figure Page 

1 Model or a fluidized bed. Particle Entrainment 15 

decreases exponentially uith increasing freeboard 
height . 

2 Side and top vieus of sampling apparatus. Z8 

3 Perspective vieu of sampling apparatus. 29 

4 Top and side vieu of sample trap. 31 

5 Design of closure paddles. 32 

S Torque output curve for rotary solenoid. 34 

7 Schematic of solenoid pouer supply. 35 

8 Schematic diagram for vacuum system. 38 

9 Heat exchanger tube design shouing the four (4) 41 

rous of 22 tubes. 

10 Heat exchanger tube design. 42 

11 Position of sample trap, bubble probe, and 44 

anemometer probe above heat exchanger tubes . 

1Z Position of bubble and Anemometer probe uith 45 

respect to the sample trap and distributor. 

13 Slock diagram of the equipment used during the 47 

sampling operations. 

14 Typical oscilloscope trace obtained during 52 

sampling operation. 

15 Maximum particle height vs particle diameter 65 

obtained during increment sensitivity analysis. 

16 Particle size distribution of bed mass used in 65 

increment sensitivity analysis. 

17 Plots of relative particle density vs freeboard 66 

height shouing the effect of varying the diameter 
interval and height interval. 



a 



18 

13 

20 

zi 

22 

23 

24 

2S 

26 

27 

28 

23 

30 

31 

32 

33 



Plot ot pressure drop across bed vs Uo. The S3 

estimate of Umf from this plot is 15,2 cm/s (0.5 
ft/s). 



Plot of particle density vs freeboard height as a 73 
function of Uo/Umf. The data uas obtained using 
the particle sampler. 

Plot of Po vs (Uo/Umf -1) shouing strong 76 

dependence of Po on Uo. 

Plot of 1/a vs Uo shouing linear dependence of 79 

1/a on Uo. 

Particle size vs mass distribution and particle 81 

number of bed material from experimental data. 

Analysis of variation of particle size vs 84 

freeboard height. 

Oscilloscope trace of bubble probe, anemometer 86 

probe, and solenoid activation at lou Uo. 

Oscilloscope trace of bubble probe, anemometer 86 

probe, and solenoid activation at lou Uo. 

Oscilloscope trace of bubble probe, anemometer 87 

probe, and solenoid activation at higher Uo. 

Oscilloscope trace of bubble probe, anemometer 87 

probe, and solenoid activation at higher Uo. 

Typical oscilloscope trace during sampling 93 

procedure shouing time before actuation. 

Relative particle nimber distributions of bed 100 

material by sieve and image analyzer analysis. 

Maximum particle height .vs particle diameter for 100 
baseline conditions. 

Relative particle number vs particle diameter for 102 
baseline conditions. Freeboard height of 4 cm. 

Relative particle number vs particle diameter for 102 
baseline conditions. Freecoard height of 8 cm. 



Reiative particle number vs particle diameter for 103 
baseline conditions. Freeboard height of 12 cm. 



7 



34 

35 

3S 

37 

38 

39 

40 

41 

4Z 

43 

44 

45 

46 

47 

48 

49 

50 



Relative particle number vs particle diameter for 
baseline conditions. Freeboard height of 18 cm. 



Relative particle number vs particle diameter for 
baseline conditions. Freeboard height of 22 cm. 



Relative particle number vs particle diameter for 
baseline conditions. Freeboard height of 31 cm. 



Particle density/unit volume vs freeboard height 
for baseline conditions. 



In particle density/unit volume vs freeboard 
height for baseline conditions. 



Maximum particle height vs particle diameter as a 
function of Uo. 



Relative particle number vs particle 
function of Uo. Freeboard height of 



diameter as a 
4 cm. 



Relative particle number vs particle 
function of Uo. Freeboard height of 



diameter as a 
8 cm. 



Relative particle number vs particle 
function of Uo. Freeboard height of 



diameter as a 
12 cm. 



Relative particle number vs particle 
function of Uo. Freeboard height of 



diameter as a 
18 cm. 



Relative particle number vs particle 
function of Uo. Freeboard height of 



diameter as a 
22 cm. 



Relative particle number vs particle 
function of Uo. Freeboard height of 



diameter as a 
31 cm. 



Particle density/unit volume vs freeboard height 
as a function of Uo. 



In Particle density/unit volume vs freeboard 
height as a function of Uo. 



Maximum particle height vs particle diameter as a 
function of Upo. 



Relative 

function 



particle number vs 
of Upo- Freeboard 



particle diameter 
height of 4 cm. 



as a 



Relative 
funct i on 



particle number vs particle diameter 
of Upo. Freeboard height of 8 cm. 



as 



a 



103 

104 

104 

105 
105 
109 

109 

110 
110 
111 
111 
112 
114 
114 
116 
117 
117 



8 



51 

52 

53 

54 

55 

55 

57 

53 

59 

60 

61 

62 

63 

64 

65 

56 



118 

118 

US 

119 

120 

120 

122 

124 

124 

125 

125 

126 

126 

127 

127 

123 



Relative particle number vs particle diameter as a 
function of Upo. Freeboard height of 12 cm. 

Relative particle number vs particle diameter as a 
function of Upo. Freeboard height of 18 cm. 

Relative particle number vs particle diameter as a 
function of Upo. Freeboard height of 22 cm. 

Relative particle number vs particle diameter as a 
function of Upo. Freeboard height of 31 cm. 

Particle density/unit volume vs freeboard height 
as a function of Upo. 

In particle density/unit volume vs freeboard 
height as a function of Upo. 

Maximum particle height vs particle diameter as a 
function of Uj . 

Relative particle number vs particle diameter as a 
function of Uj . Freeboard height of 4 cm. 

Relative particle number vs particle diameter as a 
function of Uj . Freeboard height of 8 cm. 

Relative particle number vs particle diameter as a 
function of Uj . Freeboard height of 12 cm. 

Relative particle number vs particle diameter as a 
function of Uj . Freeboard height of 18 cm. 

Relative particle number vs particle diameter as a 
function of Uj . Freeboard height of 22 cm. 

Relative particle number vs particle diameter as a 
function of Uj . Freeboard height of 31 cm. 



Particle density/unit volume vs freeboard height 
as a function of Uj - 



In Particle density/unit volume vs freeboard 
height as a function of Uj . 

Maximum particle height vs particle diameter as a 
funct i on of t j . 



3 



67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 



130 

130 

131 

131 

132 

132 

133 

133 

135 

136 

136 

137 

137 

138 

138 

139 



Relative particle number vs particle diameter as a 
function of t j . Freeboard height of 4 cm. 

Relative particle number vs particle diameter as a 
function of t j . Freeboard height of 8 cm. 

Relative particle number vs particle diameter as a 
function of tj . Freeboard height of 12 cm. 

Relative particle number vs particle diameter as a 
function of t j . Freeboard height of 18 cm. 

Relative particle number vs particle diameter as a 
function of t j . Freeboard height of 22 cm. 

Relative particle number vs particle diameter as a 
function of t j . Freeboard height of 31 cm. 



Particle density/unit volume vs freeboard height 
as a function of t j . 



In Particle density/unit volume vs freeboard 
height as a function of t j . 

Relative particle number vs particle diameter for 
bed mass material. 



Relative particle number vs particle diameter as a 
function of bed mass. Freeboard height of 4 cm. 

Relative particle number vs particle diameter as a 
function of bed mass. Freeboard height of 8 cm. 



Relative particle number vs particle diameter as a 
function of bed mass. Freeboard height of 12 cm. 



Relative particle number vs particle diameter as a 
function of bed mass. Freeboard height of 18 cm. 



Relative particle number vs particle diameter as a 
function of bed mass. Freeboard height of 22 cm. 



Relative particle number vs particle diameter as a 
function of bed mass. Freeboard height of 31 cm. 



Particle density/unit volume vs freeboard height 
as a function of bed mass. The In plot is also 
shoun as Fig. 82 b. 



10 



83 



141 



Particle density/unit volume vs freeboard height 
for baseline parameters. 

84 In Particle density/unit volume vs freeboard 141 

height for baseline parameters. 

85 Comparison of slopes for the particle density 143 

distributions above the bed as derived from 
experimental data and computer model output. 

A-l Constuction of paddles with aluminum interface 152 

cylinder shoun. 

A-Z Diagram for moment of inertia calculation used 154 

for cylinder about Z axis. 

A-3 Diagram for moment of inertia calculation used 154 

for rectangular prisim about X axis. 

C-l Diagram of closure time determination set up. 159 

C“Z Oscilloscope trace of paddle eclipsing electric 161 

eye. 

D-l Torque output of rotary solenoid showing triangle 165 

approximat ion and spring constant determination. 

M-l Calibration of anemometer probe. Oscilloscope 249 

voltage vs air velocity. 



il 



LIST OF TABLES 



Taol e Page 

1 List of equipment used during particle sampling 48 

operations . 

2 List of parameters used to check computer 62 

calculations against closed form solution, 

3 Listing of input and resulting maximum particle 64 

heights uith time to maximum height. These 

values uere used during the increment sensitivity 
tests . 

4 List of experimental data showing sample 71 

averages, standard deviations, heights, and 
velocity conditions measured. Density values are 
calculated by dividing the average sample weight 

by the sample trap volume. 

5 Results of linear, regression analysis for 75 

particle loading density (grams/cm ) vs height 

above the bed surface (cm). 

6 Comparison of least square fit relations for Po 77 

as a function of Uo and Umf . 

7 Comparison of least square fit relations for 1/a 78 

as a function of Uo. 

8 Statistical values for particle number 83 

distribution as a function of freeboard height. 

A complete listing of the data is given in 
Appendix J. 

9 Average and standard deviation of jet velocity 89 

determined from oscilloscope traces in Appendix 

K. 

10 List of transit times for particles traveling 91 

from the bed surface to the center of the trap. 

The time prior to actuation of the sample trap is 
also shown. 



1Z 



11 



35 



List of samples and their correlation parameters. 

The resulting average f or 0 indicates that no 
correlation can be made from the data obtained to 
indicate by sample weight whether or not any bed 
activity occured below the sample trap. 

12 Effect of Uo on the slope of the particle density 113 
distribution as a function of height for the 
distributions shown in Fig. 47. 

13 Effect of Upo on the slope of the particle density 121 
distribution as a function of height for the 
distributions shown in Fig. 56. 

14 Effect of Uj on the slope of the particle density 123 

distribution as a function of height for the 
distributions shown in Fig. 65. 

15 Effect of tj on the slope of the particle density 134 

distribution as a function of height for the 
distributions shown in Fig. 74. 

15' Baseline parameters used in computer model. 140 

17 Comparison of slopes for the particle density 142 

distribution above the bed as derived from the 
experimental data and the computer model. 

A. 1 Listing of paddle components and parameters. 151 

A. 2 Calculated values for the moment of inertia of 155 

paddles and paddle components. 

8.1 Results of vacuum sample removal test. 157 

C. l Sample trap closure data. 160 

D. l Summary of dynamic analysis results. 168 

E. l Average particle size distribution of bed 170 

material in grams and percentage of total weight 
using sieve analysis. 

H. 1 List of components for sampling apparatus. 179 

H.2 Listing of comoonents for solenoid power supply. 180 

H.3 Listing of components for vacuum system. 181 

distribut ions shown in Fig. 47. 



13 



CHAPTER I 



INTRODUCTION 

Fluidized beds have been used in industry for nany years. 
They have been used to mix and dry particulate materials and are 
the principle process in catalytic cracking plants. In the past 
decade or so, the use of fluidized bed combustors for power 
generation has become a source of major interest. Prototype coal 
burning beds have already been built uhich are comparable to 
existing coal plants. Fluidized bed combustors have the added 
benifit of lou NOx , SO 2 an d hydrocarbon emissions and the 
flexibility of being able to burn a wide range of fuels ranging 
from refuse and high sulfur content coal to high grade fuels. 

A fluidized bed [Fig 1] is composed of a distributor through 
uhich an air flow is introduced through thousands of small 
orfices. This air then passes through the dense zone of the bed 
uhich is comprised of a mass of particles. The air velocity 
through the dense zone is maintained above the minimum 
fluidization velocity < Umf ) during normal operations. At 
velocities equal to or greater than Umf, the frictional force 
(Drag) of the air flowing past a particle is equal to the weight 
of the particle. Under these conditions, the particle mass 
behaves very much like a fluid. It will maintain a horizontal 
surface if the container is tilted, flow out of holes in the 
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Fig. 1 Model of a fluidized bed. Particle 

entramnent decreases exponentially 
uith increasing freeboard height. 
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container, and has a pressure drop across any section of the bed 
apprxamaiely equal to the weight of air and particles in the 
section Ell. At velocities above Unf , the excess air will 
coalesce and forn small voids or bubbles as it procedes towards 
the surface of the fluidized bed. 

As bubbles break at the surface of the bed, the solid 
particles are thrown up above the bed surface and are entrained by 
the upward flowing gas stream. This zone above the bed surface is 
the freeboard zone. In the freeboard, some particles are carried 
by the gas flow far above the bed surface and are removed from the 
fluidized bed (elutriated) , while the remainder fall back to the 
bed. In general the amount of bed solids suspended in the 
freeboard (entrainment) decreases exponentially along the freeboard 
height. This distribution is similar to that of the 
Maxwel 1 -Boltzmann distribution for the case of discrete energy 
states as it applies to the Law of Atmospheres (21. 



N( z ) = No exp<-mgz/kT> 



( 1 ) 



Research in the area of entrainment by Lewis et al . [31, Zenz 
and Ueil E4] and others has resulted in the following correlation 
for entrainment as a function of gas velocity and freeboard height 
for small-particle beds [13. 



1G 



F 

At Uo 



(Z) 



Z 

= B expf - C ( b/Uo ) + a H ] > 



uhere- 



F = Entrainment 
At = Area of bed 

Uo = Superficial velocity through bed 
B = Particle dependent constant 
b = Particle dependent constant 
a = Particle dependent constant 
H = Height of freeboard 



At some point above the bed surface, the quantity of 
entrained particles becomes constant. At this point, the free 
fall velocity of the remaining particles is equal to or less than 
the uniform superficial operating velocity. The height at which 
the entrainment becomes constant is called the transport 
disengagement height (TDH) C13. 

The particles that are thrown above the bed are affected in 
the freeboard region by hydrodynamic parameters such as= bubble 
size, bubble frequency, fluidizing velocity, height above the bed, 
particle size, particle density, column diameter [S3, and 
baffles CB1. The intermittent high velocity bursts of gas which 
occurs when a bubble bursts, imposes a flucuating and highly 
irregular time dependent velocity profile over the cross section 
of the bed surface. At successively higher levels above the bed 
surface, this velocity profile becomes more and more uniform until 
at the TDH, the flow is at the uniform superficial operating 
vel oci t y ( Uo > [43. 



17 



Until recently, little attention was paid to the 
unaer standi ng of the freeboard reactions for large particle beds. 
However, due to recent work in fluidized bed combustors, the 
extent of particle and fines loading in the freeboard has been 
shown to significantly affect the S&2 absorbtion, NOx reduction, CO 

i 

emission. In general, the last 5 to 10X of the combustibles will 
burn in the freeboard. It was shown that the fine sorbent 
particles entrained into the freeboard will enhance sulphur 
capture and that entrained char particles will react with NOx and 
reduce its emission C 7 1 . Fines reinjection has been shown to 
significantly increases the fine particle concentration in the bed 
and in the freeboard with the consequence of further enhancing 
char oxidation. However, this can result in overheating .in the 
freeboard region and excess SO 2 and NOx emission C71 . The end 
result is that the potential for fluidized bed combustion power 
plants is enhanced by their ability to burn high sulphur content 
fuels and maintain low SO 2 emissions. Further research in the 
area directly above the bed surface is required to properly model 
the reactions occuring within the freeboard. 

Extensive studies on entrainment rate and elutriation have 
been made with numerous corr el at i ons , some of which are proposed 
in Cl, 3, 4, 5, 6, 7, 8]. However, most of the reported work on 
entrainment from fluidized beds has been carried out with either 
a closely sized fraction of particles or a mixture of two such 
fractions. Virtually all of this work has been conducted on bench 
scale or catalytic cracking fluidized beds. The results of these 
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studies have been shoun to produce much louer entrainment rates 
than full sized beds or tend to operate in the slugging condition. 
Therefore, the entrainment rates and transport disengagement 
heights (TDH) for fluidized beds are generally estimated from 
empirical or semiempir ical correlations obtained from this data 
and most of them shou extreme discrepancies between different 
experimental results. Extrapolation of these empirical 
correlations usually leads to strange results [51 with 
discrepancies which can vary by two orders of magnitude. 

~ The lack of good correlations stems mainly from the 
difficulty in obtaining accurate entrainment rate data. Most of 
the data is based on pressure measurements at incremental heights 
in fluidized cracking plant type beds [113. The effect of wall 
loading by particles and the actual relation between pressure and 
particle concentration is considered to be major sources of error 
when using this method with large particles. As a result, none of 
the correlations are widely accepted as giving accurate 
predictions [ 81 . 

A complete model of the entrainment process from fluidized 
beds must take into account all the mechanisims involved within 
the process. The arrival of bubbles at the bed surface, ejection 
of particles from the dense-phased bed into the freeboard region 
as the bubbles erupt, part icl e-part icle interactions, and the 
trajectories of ejected particles are all important [83. Much 
work has been done concerning bubble growth, velocity, volume, 
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etc. and their behavior is fairly uell understood. 



The mechanism of solids ejection at the bubbling bed surface 
is still not uell understood. The origin of ejected particles is 
reported to be primarily due to tuo sources. The particles uhich 
have been lifted by the bubble uake and thrown upwards following 
the bubble burst at the surface is the -first source. This theory 
is supported by work done by George and Grace (81 who performed 
experiments which concluded that the vast majority of the ejected 
particles did not originate from the surface layers but from 
bubble uake pick up. Uork done by._Page and Harrison CS1 also 
appears to agree with this* The second theory suggests that the 
ejected particles originate at the nose of the bursting bubbles 
and are thrown outward when the bubble breaks. Research by Rowe 
and Partridge [8] and Glicksman et al [12] have shown this second 
mechanisim as being the dominate particle ejection source and thus 
supporting this second theory. Their work has also shown that 
under the conditions in which 2 bubbles coalesce just below the 
surface of the bed, the jet of gas produced can result in a 
significant amount of particles being ejected from the wake of 
the first bubble. 

The effect of mu 1 1 i part i cl e interactions have been for the 
most part ignored except by Peters and Prybylouski [133. The 
motion of any individual particle is influenced by the presence 
of other particles, i.e., through direct part i cl e-par 1 1 cl e 
interactions and deviations in the fluid drag force. The major 
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drawback of their work is that the paper compares their theory 
with only a single set of experimental results C 31 - 

Studies to model the trajectories of particles in the 
freeboard have been conducted several times. The uork of Ualsh 
- et al^L7], George and Grace E8] t and Peters and Prybylouski [131 
name just a feu of_the latest efforts. All of these studies 
relied upon experimental data to develope their theories* 

Houever , to check the accuracy of their theories, more 
experimental data is required. 

As of yet, none of the entrainment models available can be 
incorporated into fluidized bed combustion models uith sufficient 
accuracy to uarrent their use. This is due to a lack of 
experimental information on entrainment rate as a function of the 
complete fluidization parameters of the bed to test the models 
uith. As a result, the purpose of this study has been to obtain 
particle density distributions above a cold atmospheric fluidized 
bed containing a continuous particle size distribution [Appendix 

El . 
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CHAPTER XX 



PARTICLE SAMPLING APPARATUS 



Design Alternatives 

There are many methods available for determining particle 
distributions in fluid flous. The more commonly used methods are*’ 

1) Catching mechanisms 

Z) Trapping mechanisms 

3) Radiation attenuation measurements 

4) Optical measurements 

5) Capacitance and Inductance measurements 

Catching mechanisms are passive devices. That is, particles 
are captured merely by the presence of the catching mechanism in 
the fluid flow containg the particles to be sampled. The data 
obtained using this method is position dependent and produces 
average values for the particle flux loadings. These catching 
mechamsims are also limited in that they can only catch particles 
with particle fluxes traveling in a single direction. The device 
used by Ualsh et al [10] only captured falling particles uhile the 
device used by George and Grace [8] required the uouard moving 
particles to deflect off of a baffle surface and fall into a 
collecting trough. 
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Trapping mechanisims, unlike catching mechanisms, are active 



particle samplers. Their operation involves the trapping and 
isolation of a finite volume of the fluid flow at a specific 
period in tine. This sampling technique produces time dependent 
as well as position dependent data. This uill allou correlations 
between bubble eruption and particle density to be made using 
multiple bubble conditions rather than single bubble capture. As 
the number of random samples taken by this method increases, the 
average value of this data uill approach that of the catching 
mechanism. Trapping mechanisims also capture particle fluxes 
traveling in multiple directions- This ability reduces the error 
inherent in measuring only the downward or only the upward 
particle flux. The apparatus used in this paper is a trapping 
mechanism. 

Attenuation of nuclear particles from a radioactve source can 
be used to give average particle density distributions across a 
suspension. However, this method is not adaptable to density 
determi nat i ons at a point. An average time dependent density 
determination can be achieved with this method. Another draw-back 
of this method is the radiation hazards involved with the use of 
nuclear particles. 

Optical density determinations consists of two seperate 
methods. The first method uses a very small light beam which is 
eclipsed by the transition of a particle through it. A related 
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method uses the absorbtion and scattering of a somewhat larger 
light beam to correlate the change in light intensity with 
particle density. This method has been used frequently in the 
study of aerosols but requires complicated and intricate equipment 
C141. The second method involves high speed photographs of a 
small volume of space. This method cannot be used uhen the 
particle density. is^so large that multiple particles eclipse each 
other frequently enough to produce unacceptable error. This is 
the case uhen the probe height above the bed is less than 7-15 cm 
(3-S in). 

~ ' Measurements at a point can also be made by inserting either 

a toroidial inductor or a parallel plate capacitor in the flow. 

The presence of the particles changes the permeability and thus 
the inductance of the inductor, or the dielectric strength and 
thus the capacitance of the capacitor. The drau-backs of these 
methods involves the unknoun effects of particle velocity and 
external particles on the inductor and charge transfer to 
particles from the capacitor (141. 

Apparatus Requirements 

The goal of this study uas to determine the density 
distribution of particles above an atmospheric fluidized bed uith 
particle velocities of up to 10 meters per second. The particle 
size distribution of material ejected from bubbles is required for 
particle trajectory calculations. A correlation betueen the 
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average density and the density present immediately after a bubble 
bursts fron the bed surface uas also of interest. These 
requirements dictated that the method used for measuring densities 
have the following capabilities^ 

1) Measure densities with good spatial resolution. 

Z) Measure densities at specific moments in time. 

3) Obtain particle size information. 

4) Operate under extremely dirty conditions. 

5) Easy sample removal from bed. 

G) Remote operation of sampler. 

The radiation attenuation and inductance/capaci tance methods 
can not determine particle sizes- Therefore, these methods were 
no longer considered as possible measurement alternatives. 

Because the optical methods are either not reliable at small 
heights above the bed or their use is too complex, they uere not 
used. Catching devices, although simple to use, do not have the 
ability to measure data at specific points in time and determine 
particle density loading in space. As a result, the determination 
to use a trapping mechanism as the method of measuring particle 
densities uas made. 
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Apparatus Design 



General Design Criteria 

The following criteria was used in determining the design of 
the trapping device. 

1) The closure tine of the trap was chosen to be equal 
to the tine for a particle with a velocity of 10 
n/s to transit 1/10 the length of the sample 

~ container. This velocity is considered to be the 
upper limit of the particle velocity distribution 
present in the test bed, based on the work of 
George and Grace C8], 

2) The apparatus must be capable of frequent sampling 
without requiring access to the sampling device 

i tsel f . 

3) The samples trapped, must be easily accessable 
from outside the fluidized bed without interupting 
the bed conditions. 

4) The apparatus must be able to operate in the high 
particle flux environment of the fluidized bed. 
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5) The actuation of the trapping device must be able 



to be accurately determined to allou correlation 
with other time resolved measurements. 



General Design 

The apparatus is shoun in Figs. 2 and 3. A description and 
list of all components is given in Appendix H. The sample 
container is mounted on an extension arm to minimize the 
^disturbance to the air flow around the sample trap caused by the 
rest of the mechanism. The sample container is closed using two 
(2) paddle arms, one above and the other belou. These paddle arms 
are attached to aluminum interfaces which are used to connect them 
to a rotary solenoid. The solenoid is used to suing the paddles 
over the sample trap and shut it. Not shoun in these figures are 
the power supply for the solenoid, the vacuum system used to 
remove the particles from the sample trap and the water-proof 
nylon shell used to keep the particles from interfering with the 
operation of the solenoid. All of these systems are described in 
greater detail in the following sections. 

Sample Container 

To ensure that the sampling device had minimal effect on the 

fluid flow, the cross sectional area presented to the flow had to 

be minimized. This constraint required that the sample container 
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Fig. 3 Perspective view of oampltnQ apparatus 



be seperated fron the rest of the apparatus. This uas benificial 
in the final design because it helped to reduce the apparatus 
closure tine. 

To ensure a good seal uas achieved uhen the trap uas shut, 
felt uas used as a gasket between the sample container and the 
closure paddles. Fig. 4 shows thehf i nal design used for the 
sampling mechanism. It is constructed of 1/1S inch aluminum with 
1/4 inch square stock used for the frame and mounting structure. 
Epoxy is used to seal the sides of the container. 

Closure Paddles 

After several iterations on paddle designs, it uas determined 
that the paddle construction which offered the greatest stiffness 
for the least weight uas a composite laminate. The paddle, shown 
in Fig. 5, is made using a 0.4 in thick foam core with 1/32 inch 
thick Basswood laminations on the exterior. Harduood (Maple) end 
pieces were used to provide a noncompr essi ve connection between 
the foam paddles and the aluminum interfaces. The aluminum 
interfaces couple the solenoid shaft to the paddles. Epoxy uas 
used to join the laminate materials. 

Actuator 

A rotating mechanism utilizing a rotary solenoid uas chosen 
to shut the sample trap. A rotary solenoid was selected because 
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of the problems of maintaining low friction surfaces for sliding 
mechanisms in the presence of the particulate matter being 
sampled. This precluded the use of any sliding mechanism to shut 
the sample trap. The calculations found in Appendix A and 
Appendix C determined the size of the solenoid required to achieve 
the desired closure time. 

The 45 degree stroke solenoid was chosen to place the paddle 
arms far enough auay from the sample trap, such that when 
de-energized it prevents interference with the particle flou. 

This stroke also minimized the area which must be clear of 
— obstruct ions-to -the travel'of the paddles. The solenoid operates 
at a 1/10 duty cycle pouer rating uhen initially actuated, 
providing the torque output shown in Fig. 6. After the solenoid 
has shut the sample trap, the solenoid is operated at a lower 
power rating, providing a holding torque of 5.5 in-lbs. This 
decreased rating is necessary to prevent overheating of the 
solenoid. This assembly is encased within a nylon shell. 

Power Supply 

Fig. 7 is a schematic of the electrical system used to power 
the solenoid. Appendix H contains a list of all components used 
in the power supply. The power supply plugs directly into a 
standard 115 volt AC line source. Switch SI is used to apply 
power to the solenoid Ml. The full wave bridge rectifier assembly 
converts the AC line voltage to DC. The rectifier assembly is 
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Schematic of Solenoid Power Supply. 
Q Heference number for component 



protected by the arc suppresser 01 from the large voltage spike 



induced in the solenoid coil uhen the solenoid is actuated. 
Resistor R3 is used to reduce the current through the solenoid 
after the solenoid has been closed to prevent it from overheating,, 
To initially shut the sanple trap, full current is applied to the 
solenoid by shutting suitch SZ , uhich bypasses resistor R3. 
Resistors R1 and RZ forn a voltage divider netuork to provide a 
low voltage (5.5 volts AC) trigger source for an oscilloscope. 
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Vacuum System 



To remove the particles in the sample trap, a vacuum system 
was developed. This system is shown in Fig. 8 uith a list of the 
components given in Appendix H. A vacuum is produced by allowing 
air from a 100 psi air source to flow through valve VI into the 
venturi eductor PI. The vacuum places a suction on the sample 
container CZ via a fine mesh screen. The purpose of this screen 
is to prevent particles from escaping the sample container. The 
suction is applied to the sample trap Cl through 1/4 inch polyflow 
tubing. It is through this tube that the particles are removed 
from the sample trap and collected in the sample container. An 
equalization and agitaion line is connected to the opposite side 
of the sample trap. This line serves two purposes. First, it 
ensures that the vacuum system does not pull in particles from 
outside of the sample trap. Second, it allows a flow of air to be 
introduced which stirs up the particles trapped inside. • This 
helps push them into the suction line and reduce the remaining 
particles to a minimum. 

Apparatus Testing 

Two tests were run to determine the effectiveness of the 
system. The first test determined the closing tine of the sample 
trap. The second evaluated the error from the loss of particles 
which were left in the sample trap by the vacuum system. 
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Schematic DiaQram For Vacuum System. 
O Reference number for component 



The procedure and method used to determine closure time for 



the sample trap is given in Appendix C. From the results of these 
tests given in Appendix C, the closure time uas determined to have 
an average value of 1.44 milliseconds. This is equal to a 10 m/s 
particle traveling 1.44 cm or approximately IS Z of the sample 
trap length. The average particle velocity uill be less than Z 
m/s and uill therefore introduce an average error of less than 3 Z. 

Using the sample trap closure time, a dynamic analysis 
determined that the total time from initially applying power to 
the solenoid until it shut the sample trap is 4Z.S milliseconds. 
These calculations are given in Appendix D. This actuation time 
is important for determining uhether a specific bubble’s debris 
uas within the vicinity of the trap at the time of closure. 

The procedure and results for the testing of the sample 
removal vacuum system are given in Appendix 3. The average amount 
of particles lost by the vacuum system uas determined to be 0.5Z Z 
of the initial sample placed in the trap. The maximum error uas 
0.33 Z. The error from the vacuum system is therefore considered 
to have an insignificant effect on the data obtained. 
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CHAPTER III 



EXPERIMENTAL PROCEDURE 
Fluidized Bed Conf i cur ati on 

The M.I.T. atmospheric fluidized bed, in which the sampling 
device was used, is a model of the 20 MU atmospheric fluidized bed 
combustor prototype, jointly sponsored by the Tennessee Valley 
ftjihority and the Electric Pouer Research Institute. The 
fluidized bed model is described in Lord et al (151 and Jones et 
al £163. While using the sampling apparatus, a different heat 
exchanger tube bundle conf igurat ion uas used than is described in 
Jones et al [163. The heat exchanger conf igurat i on used is shoun 
in Figs. 9 and 10. 

The heat exchanger used during this work is made of 1.26 cm 
(0.5 in) O.D. tubing arranged in 4 rous of 22 pipes each. The 
tubes are aligned as shoun in Fig. 10. Each pipe is spaced uith 
a vertical center to center distance of 5.08 cm (2 in) and a 
horizontal center to center distance of 3.91 cm (1.5 in). The 
distance from the distributor to the center of the upper most tube 
is 27.52 cm (10.875 in). A 5.08 cm (2 in) spacing separates the 
front and back walls of the fluidized bed from the end tubes of 
the bundel. The cross sectional area of the bed is 1.079 sq m 
(11.61 Sq Ft). The particulate material used in the bed is a 
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Fig. 9 Heat Exchanger Tube Design Showing The 

Four (4) Rous of ZZ Tubes. 

All Dimensions in cn. 
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mixture of steel grit abrasive having a specific gravity of 8.1. 
Appendix E lists the size distribution of the steel grit used 
during the sampling operation. The bed uas operated without 
recycling the fines captured in the cyclones. The static bed 
height of the material uas 22.54 cm (8.875 in) throughout the data 
collection period. 

Equipment Set-ug 

Fig. 11 shows the placement of the sampling device inside 
the fluidized bed. The sample trap uas positioned such that it 
uas directly over a spacing between tubes 122.36 cm (9 in) from 
the center of the front-most tube] and 41.91 cm (16.5 in) from the 
left uall. The height of the sample trap above the distributor 
uas varied during the sampling process as discussed in the section 
on sampling procedure. 

Fig. 12 shous the placement of the bubble probe and the 
anemometer probe uith respect to the sample trap. The bubble 
probe uas placed directly below the sample trap and 26.67 cm (10.5 
in) above the distributor. The probe extension was placed at an 
angle so as not to interfere uith the sample trap operation. To 
protect the anemometer wire from particles impacting it, a special 
shield consisting of S 320 mesh screen and an aluminum frame uas 
placed around it. The anemometer probe uas attached to the bubble 
probe extension uith the entrance to the anemometer probe 29.21 cm 
(11.5 in) above the distributor. This placed the entrance to the 
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anemometer probe 2.54 cm (1.0 in) auay from the center of the 
sample trap and 0.S4 cm (0.25 in) outside the area directly belou 
the sample trap. As a result, the probe has a minimum effect on 
the air flow from the bed to the sample trap, but will only 
measure the gas velocity at the edge of the sample trap perimeter. 
The positioning of the two probes above the distributor remained 
constant throughout the sampling evolution. 

Fig. 13 is a block diagram showing the equipment used during 
the sampling operations and their 'interconnections . Table 1 is a 
listing of the equipment used. The oscilloscope time base uas set 
for MANUAL TRIGGER, SINGLE SUEEP mode and a sweep time of 50 
ms/div. The channels of the dual trace amplifier were set at 5 
volts per division for the bubble probe and 2 volts per division 
for the sample trap inputs. On the differential, amplifier, one 
channel uas not used and the second channel uas set at 1 volt per 
division for the anemometer probe input. 
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Fig. 13 Block Diagram of the Equipment Used 
During the Sampling Operations. 
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TABLE 1 



Osci 1 loscope 


TEKTRONIX 5111 Storage Oscilloscope 
TEKTRONIX 5A18N Oual Trace Ampl 
TEKTRONIX 5A21N Differential Ampl 
TEKTRONIX 5B10N Time Base 


Anemometer 


Thermal Systems Inc. 

1051-Z Monitor and Pouer Supply 
10S4-A Linearized Anemometer 
Anemometer uire u/ £323 screen guard 


Optical Signal 


Detector 



Optical Signal Detector Pouer Supply 
Optical Source and Pouer Supply 
Oscilloscope Camera 

Listing of equipment used during particle sampling 
operations . 
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To determine the fluidization conditions uithin the bed, a 
set of nanometers were used. These manometers measured pressures 
uithin the bed, at heights from 4.13 cm (1.62 in) to 37.15 cm 
(14.52 in) above the distributor, in 2.54 cm (1 in) increments. 

The pressure data correspond! ng to each trap position and bed 
velocity is listed in Appendix I. 

To determine the gas flou conditions uithin the bed, an 
orfice flou meter uith ID - 1/2D taps uas located upstream of the 
distributor. The computer program listed in Appendix G uas used 
to convert the pressure tap data to mean air velocities uithin the 
bed. 

Sampling Procedure 

Data uas collected for four (4) mean bed velocities at six 
(6) different sample heights. The sample trap uas placed at a 
given sampling height (measured from the distributor to the bottom 
of the trap), and ten (10) samples uere collected at each desired 
velocity. The trap position uas then changed to a neu height. 

During certain sampling conditions, those unich involved lou 
sampling heights uith the higher air velocities, the paddle arms 
uould occasionally impact the sides of the sample trap and not 
close the sample trap completely. It is assumed that this occured 
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uhen a large bubble erupted directly under the paddle arm and 
deflected the paddle arm into the side of the trap. Whenever this 
occured, the trap was de-energized and the closure cycle repeated. 

Each time the height or velocity was changed, a complete 
sampling cycle uas conducted and this sample discarded. This uas to 
prevent any accumulation of particles (in the entrance to the 
vacuum or purge lines on the sample trap) from being added to the 
first sample at the neu height or velocity. 

For each set of data at a given height and velocity the 
following information uas recorded- 

1) Fluidized bed height above the distributor 

determined visually and by pressure measurements. 

Z) Pressure upstream of orfice plate (PI) 

3) Pressure difference across orfice plate (delta P) 

4) Air temperature in bed 

5) Pressure distribution in bed 

G) Height of sample trap above distributor 

The following procedure was used during sampling^ 

For each sample to be taken within a data set. 

A. Initial conditions 

1) Vacuum air supply OFF 

2 ) Sweep trigger on oscilloscope RESET 
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3) O-scope memory 



ON-CLEAR 



4) Solenoid pouer supply switch Si 



OFF 



5) Solenoid trigger switch SZ 



ON 



B. Sampling Procedure 



1) Trigger oscilloscope sweep and wait until sweep 
is at the center of the CRT, 

Z) Close the Solenoid power supply switch Si. Uhen 
SI is shut, the oscilloscope will show an 
additional trace. This third trace is used to 
determine the closure time relative to the 
presence of gas jets and bubble eruptions. An 
example of a typical oscilloscope trace is shown 
in Fig. 14. 

3) Open the solenoid trigger switch SZ. This 
reduces the current to the solenoid to prevent 
overheating. The maximum allowed time to let SZ 
remain closed is five (5) seconds. 

4) Turn on the air supply to the vacuum system and 
leave cn one (1) minute. The exhaust air from 
the venturi on the vacuum system must be directed 
into the purge line in an oscillatory manner. 

This will agitate the particles within the trap 
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Fig. 14 Typical Oscilloscope Trace Obtained 
During Sanpling Operation. 
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so as to move then into the vacuum line. 



5) Turn off the vacuun air supply. 



S) Turn off the solenoid power supply switch Si. 



7) Renove sanple from sanple container and place in 
specinen bottle. 



8) Photograph trace on oscilloscope. 



9) Repeat fron A. 



Sanple Analysis 

Sanple Ueight Deterni nation 

Each sanple obtained was weighed, using a Torsion Balance Co. 
TORBAL scale, to an accuracy of 0.01 grans. The average value and 
standard deviation was then determined for each set (specific 
height and velocity) of samples. The weight of particles in a 
completely filled sanple trap was also determined for void 
deterni nai i on. These results are listed in Appendix I. 
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Determi nation of Particle Size Distribution 



To determine the particle size distribution uhich occurs at 
various heights for a specified gas velocity, three (3) samples 
from each set of data at a specified height and velocity, and from 
the bulk bed material were analyzed using a Zeiss Videoplan 2 
Image analyzer. The software used was the “Image Analysis System 
MQP-Videoplan ,, distributed by KONTRON Electronics Group. 

For each sample selected, a microscope slide was prepared 
using double sided adhesive tape on which a portion of the 
selected sample was placed. To ensure a sharp contrast was 
achieved, each slide was backed with white paper. The sample was 
then placed under a Zeiss microscope to which the image analyzer 
was connected via a vidiocon tube. The magnification used was 
50X uhich provided an average view of about 8 particles at a 
time. The analyzer was then used to analyze the partices present 
on each slide of interest. The mode in which the image analyzer 
was used was the equivalent diameter mode. This mode determined 
the cross sectional area of each particle sampled and calculated 
the diameter of a circle with the same area. Uhen a data set was 
completed, the data was analyzed for particle size distribution. 
The resulting output (Appendix J) consists of a particle count vs 
particle size histogram, a cumulative frequency plot and a 
classification data list. A gaussian distribution fit for the 
data is also plotted on the histogram and cumulative frequency 
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plots. The X axis of the plots are normalized uith a range of 
zero (0) to four (4). To determine the actual diameter of the 
particle in microns for a given normalized value, the normalized 
value must be multiplied by the conversion factor 248.887. 
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CHAPTER IV 



COMPUTER MODEL 



Introdjcii on 

This chapter descibes the theory, logic flow and testing of a 
particle trajectory computer model to predict the particle loading 
in the freeboard due to erupting baubles. The program, listed in 
Appendix L, is written in HP SASIC 2.0 and uas run on a Hewlett 
Packard 981S mi crocompu ter . The output from the model is discussed 
in chapter VI and compared with experimental results. An in depth 
analysis of the program logic and structure is given in Appendix L. 

Model Theory 

The model developed here, is based on calculating the 
trajectory of a single particle as it is ejected from the bed 
surface and is acted upon by gravitational and drag forces. The 
drag force is due to the difference in absolute particle and air 
velocities. The air velocity is a combination of the initial jet 
velocity produced when a bubble bursts, and the superficial bed 
velocity. To ensure that the particle drag is calculated 
accurately, the following drag coefficient correlation given by 
Uhite [171 was used. 
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Cd = Z4 + 
Re 



+ 0.4 



< 3 ) 



6 

i+ {rP 



where: 

Cd = Drag coefficient for sphere 
Re = Reynolds number 



Eqn 3 is valid over the range 0<Re<10 • To calculate the 
particles position and velocity, the computer uses a forward 
difference method- Using Newton’s Law (ZF-ma) the acceleration of 
the particle due to gravity and drag is determined- Inserting 
this acceleration into Eqn 4, the particles new velocity is 
determined. 



V = Vo + a t 



(4) 



where : 

V 

Vo 

a 

t 



Particles new velocity 
Particles present velocity 
Acceleration of particle 
Time increment of calculation 



To determine 
calculated in Eqn 



the particles 
4 is inserted 



new position, 
into Eqn 5. 



the velocity 
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H = Ho + (Vo ± V) t 
Z 



( 5 ) 



where* 

H = Particles neu height 
Ho = Particles present height 
Vo = Particles present velocity 
V = Particles neu velocity 
t - Tine increment of calculation 



Using Eqns 3,4, and 5, the trajectory of the particle is calculated 
fron the tine it initially leaves the bed until the tine that it 
returns to the bed. 

These calculations are repeated over a range of particle 
dianeters fron 80 to 570 nicrons. By determining the residence 
tine of each particle uithin a specified height increnent (A H) 
above the surface of the bed, a particle density distribution 
above the bed is deternined. The height increnent (A H) used in 
the progran is 2 cm. At the end of each tine step uhen the height 
calculation ( Eqn 5) is completed, the counter representing the 
particular Z cn height increnent uhich the particle is in, is 
incremented by one. Each particle size has its own set of 
counters to allow individual particle analysis. 

The calculated density distribution is then weighted with the 
particle size distribution of the bulk bed material since the 
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probability of a given size particle being present at a specified 
height is dependent upon the number of particles uithin the 
system. This is accomplished by multiplying each height counter 
of a given particle size uith the number of particles for that 
given size present in the input bed distribution. By summing the 
density values for each set of particle diameters at a given 
height over the entire freeboard of the bed, the overall particle 
density above the bed surface is determined. 

The model assumes that all the particles are ejected 
perpendicular to the surface of the bed and are initially at a 
uniform velocity. Because the model uses single particles for the 
analysis, the effects of multiparticle interactions are not 
included in the model. 

Testing of Proaran 



To evaluate the validity of the program, tuo tests uere 
conducted. The first test compared the height solution produced 
by the computer uith a closed form solution. The second test 
involved running the program uith different particle diameter and 
particle distribution height intervals to ensure that a valid 
sample size uas being used. 

Closed Form Solution 

To determine a closed form solution for particle height as a 
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function of initial particle velocity and superficial velocity Uo, 
a force balance uas used- The forces acting on a particle are 
gravitational and drag. The gravitational force, FI, is singly 
the volume of the particle multiplied by the particles density and 
the gravitational acceleration, and can be uritten as- 



3 

FI - - Pc IT 0 g 
S 



( 6 ) 



uher e : 



FI - Force due to gravity 
Pp = Density of particle 
D = Diameter of particle 
g = gravitational acceleration 



In order to get a closed form solution that did not involve 
non-linear differential equations, Stokes flow uas used for the 
closed form solution only. The computer model used the Stokes 
equation only to compare results uith the closed form solution, 
afterward, Eqn 3 uas used. Using the Stokes drag coefficient 
relation, the drag force on a particle can be determined as: 



F2 = 3 (Uo - Up) ^ D 



( 7 ) 



uher e : 



F2 = Force due to drag 
Uo = Superficial bed velocity 
Up = Velocity of particle 
j* = Absolute viscosity of air 
D = Diameter of particle 
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By inserting Eqns S and 7 into Newton’s Law (Z F=ma> and 
simplifying, the following differential equation is obtained 1 



X + Cl X = C2 (8) 

where* 

• * 

X = Acceleration of particle 

X - Velocity of particle 

Cl = 18^t 
Pp D 2, 

C2 = Cl Uo - g 

boundary conditions^ 

1 ) t=0 X=Uo 

2) t = 0 X=0 



This second order linear differential equation can be 
solved using the given boundary conditions uith the resulting 
closed form solution given as - 
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( 9 ) 



X = C3 + CZ t - C3 exp< -Cl t > 
Cl 



utiere- 



X = Height of particle 



C3 = 





JL 

Cl 



Using the initial conditions listed in table Z, the solution 
obtained using the computer model (maximum height^ 15.Z70 cm, time 
to maximum height^ 0.Z75 sec) uas identical to three decimal 
places with the solution obtained using Eqn 9. 



TABLE Z 



Superficial Velocity ( Uo ) - 


60.96 


cm/s (Z ft/s) 


Initial particle velocity ( Upo ) - 


304.8 


cm/s ( 10 ft/s) 


Particle diameter: 


Z00 


mi crons 


Particle density- 


5000 


kg/m^ 


Time increment (computer): 


0.001 


sec 


Jet velocity: 


0.0 


cm/s 



List of parameters used to check computer 
calculations against closed form solution. 
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Sample Size Sensitivity Test 



To ensure that appropriate sanple sizes u ere used to minimize 
errors die to coarse sampling intervals, tuo sensitivity tests 
u ere run. One test involved changing the particle diameter 
interval from 10 microns to 5 microns. The second test changed 
the height sampling interval (A H) from 2 cm to 5cm and then 1 cm. 
For each of the tests, the same initial conditions uere input into 
the program. Table 3 shous the resulting output from the program 
listing the initial conditions and the resulting diameter versus 
maximum height data calculated. Fig. 15 shous a plot of the 
calculated maximum height vs particle diameter data listed in 
table 2. Fig. 16 shous the bed particle distribution used for 
each of the tests. 

Fig. 17a shous the entrainment calculation using a diameter 
interval of 10 microns and a A H of 2 cm. Fig. 17b shous the 
same calculation using a diameter interval of 10 microns and a A H 
of 5 cm. The curve is not as smooth but still retains the same 
general shape. The peak of the curve shoun in Fig. 17b occurs at 
a height of about 20 cm uhereas the peak in Fig. 17a occurs at 
about 16 cm. A semi-log plot of these curves uould shou that the 
slope of the line to the right of the peak uould be larger for the 
data represented by Fig. 17b. The effect of maintaining A H at 5 
cm but decreasing the diameter interval to 5 microns is shoun in 

Fig. 17c. There is no readily detectable difference betueen Fig. 
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cn/s 

Mean 3ed Veiocity= 57.312 

Initial Particle Velocity 3 37.2312 
Peak Jet Veiocity= 503.5 

Gas Jet Duration* .02 s 



Diameter 


Maxinun 

Heicht 


Max Ht 
T ine 


Diameter 


Maxi nun 

Heiaht 


Max Ht 
T ime 


um 


cn 


seconds 


um 


cm 


seconds 


30 


55 . 1 


.292 


90 


57.0 


.290 


100 


57.2 


.297 


110 


56 .8 


.302 


123 


55.8 


.305 


130 


54.4 


. 307 


140 


52.7 


.307 


150 


50.3 


.305 


160 


48.8 


.303 


170 


46.8 


.299 


130 


44.7 


.295 


190 


42.7 


.290 


200 


40.7 


.285 


210 


38.8 


.280 


220 


37.0 


.275 


230 


35 . 3 


.270 


240 


33.7 


.264 


250 


32.2 


.259 


260 


30.8 


.254 


• 270 


29.5 


.249 


280 


28.2 


.245 


290 


27.1 


.240 


300 


25.1 


.236 


310 


25.1 


.232 


320 


24.2 


.228 


330 


23.3 


.224 


340 


22.5 


.220 


350 


21.3 


.217 


360 


21 . 1 


.213 


370 


20.4 


.210 


380 


19.3 


.207 


390 


19.2 


.204 


400 


18.7 


.201 


410 


18.2 


. 199 


420 


17.7 


.196 


430 


17.3 


.194 


440 


16.9 


.191 


450 


16.5 


. 189 


460 


16.1 


. 187 


470 


15.7 


.185 


480 


15.4 


. 133 


490 


15.1 


. 181 


500 


14.9 


.179 


510 


14.5 


. 178 


520 


14.2 


.176 


530 


14.0 


. 174 


540 


13.7 


.173 


550 


13.5 


. 171 


550 


13.2 


.170 


570 


13.3 


. 153 



Table 3 Listing of Input and Resulting Maximum 
Particle Heights with Tine to Maximum 
Height. Thes Values Uere Used During 
the Increment Sensitivity Tests. 

64 



Hulatlve Particle Number y Max Pari 1 cl a Height (cm) 



50 
40 
30 -■ 



10 T 



‘ ' ^ ■ 1G0 ' 200 300 ' 480 1 580 ! 

I 

1 Particle Diameter (microns) 

15 Maxirum Particle Height vs Particle 
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Fig. 15 Particle Size Distribution of Bed Mass Used 
in Increment Sensitivity Analysis. 
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Relative Particle Density T Relative Particle Density 







Plots of Relative Particle Density vs Freeboard 
Height shouin the effects of varying Diameter 
Interval and Height Interval . 
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Fig. 17 



17c and Fig. 17b. In Fig. 17d, the effect of changing A H to 1 
cn and dianeter spacing to 5 microns is shoun. The difference 
between Fig. 17a and Fig. 17d is barely noticable and no 
detectable change in the slope to the right of the peaks is 
present. As a result of this analysis, the program uas operated 
with a A H of Z cm and a diameter interval of 10 microns. This 
reduced the calculation time to half of that required uhen using a 
diameter spacing of 5 microns uhen the same total diameter span 
uas used and, as uas seen in Figs. 17a and 17d, the difference in 
output does not require the finer increment. 
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CHAPTER V 



EXPERIMENTAL RESULTS AND DISCUSSION 
Minimum Fluidization Velocity 

Fig. 18 is a plot of nean bed velocity ( Uo ) vs pressure drop 
through the bed (Pb). The velocity values uere deternined using 
the program in Appendix G and the pressure data listed in Appendix 
I- From this plot, the mininun fluidisation velocity ( Umf ) for 
the bed conditions used during this study is determined to be 0.15 
m/ s (0.5 Ft/sec). 

Entrainment Analysis 

Table 4 lists the averaged sample weights and their standard 
deviations for the samples (Appendix I) collected by the sampling 
apparatus. Looking at the standard deviation of the sample 
groups, the standard deviation is fairly large compared to the 
average values. However, visual observations of the fluidized bed 
in operation would suggest that a larger standard deviation would 
be expected. The short sample cycle (2 ms) and the bubble burst 
activity in the bed are the main reasons for this conclusion. The 
average standard deviation is 29 7 . of the average sample weight 
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( ja^en uo ) dojQ 



PIq. IQ Plot of Proaeura Drop flcrosa 0ad va Wo. The 
Eat ina to of Uwf fron thlo Plot io 15.2 cn/e 
<0.5 ft/fi). 



and the range extends from S3 I of the average for the sample set 
numbers 106 thru 115 to only 7 X of the average for the samples 
285 thru 295. This variation is believed to be due to the nature 
of the bubble activity uithin the bed. It would appear that when 
a bubble erupts at a time just prior to closure of the trap, a 
sample weight greater than the average would be obtained. 

However, when no bubble has been present, the material caught by 
the trap should consist mainly of particles returning to the bed 
and therefore be below the average weight caught . Table 4 also 
shows the density of each sample collected. These values were 
determined by dividing the averaged sample weights by the volume 
of the sample trap which is 43.02 cubic cm. 
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TABLE 4 



Superficial Bed Velocity (Uo): 58.1 cm/ s (1.905 ft/s) 

Non-dimensionalised velocity ( Uo/Umf ) : 3.81 



Samp! e 


Height 


Average 


Standard 


Average 


Numbers 


Above Bed 


Weight 


Deviation 


Density 




( cm ) 


( grams ) 


( grams ) 


(grams/cm ) 


2 85-235 


3.3^1 




0.51 


1.68 (-01) 


246-255 


7.6 


4.48 


0.72 


1.04 (-01) 


116-125 


12.4 


Z .41 


0.76 


5.60 (-02) 


166-175 


17.8 


1.08 


0.16 


2.51 (-02) 


76-85 


ZZ.Z 


1.3Z 


0.23 


3.07 (-02) 


176-185 


31.1 


0.31 


0.10 


7.21 (-03) 


Superficial 


Bed Velocity 


( Uo ) : 


48.3 cm/s 


(1.585 ft/s) 


Non-dimensionalised velocity ( Uo/Umf ) : 


3.17 




Sample 


Height 


Average 


Standard 


Average 


Numbers 


Above Bed 


Weight 


Deviation 


Density 




( cm ) 


( grams ) 


( grams ) 


(grams/cm ) 


276-285 


5.1 


4.35 


1.05 


1.01 (-01) 


236-245 


8.9 


1.88 


0.22 


4.37 (-02) 


96-105 


13.6 


1.36 


0.32 


3.16 (-02) 


156-165 


19.1 


0.60 


0.12 


1.39 (-02) 


56-65 


24.1 


0.3Z 


0.36 


2.14 (-02) 


136-195 


31.8 


0.18 


0.04 


4.18 (-03) 



Superficial Bed Velocity ( Uo ) : 
Non-dimensional ised velocity (Uo/Umf): 


39.5 cm/s 
Z . 59 


(1.Z97 ft/s) 


Sample 


Height Average 


Standard 


Average 


Numbers 


Above Bed Weight 

(cm) (grams) 


Deviation 
( grams ) 


Density 
(grams/cm ) 


266-275 


5.7 1.31 


0. 37 


3.05 (-0Z) 


226-235 


9.5 0.62 


0. 17 


1.44 (-0Z) 


126-135 


14.3 0.28 


0.08 


5.51 (-03) 


136-145 


19.7 0.08 


0.04 


1.86 (-03) 


86-95 


24.8 0.09 


0.05 


Z . 09 (-03) 


196-205 


32.4 0.06 


0.0Z 


1.39 (-03) 
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TABLE 4 ( cont ) 



Superficial Bed Velocity ( Uo ) : 35.4 cm/s (1.1S1 ft/s) 

Non-dimensional i sed velocity (Uo/Umf): 2.32 



Sampl e 


Height 


Average 


Standard 


Average 


Numbers 


Above Bed 


Ueight 


Deviation 


Density 




( cm ) 


( grans ) 


( grams ) 


(grams/ cm ) 


255-265 


6.4 


0.63 


0.13 


1.45 (-02) 


215-225 


10.2 


0.32 


0.04 


7.44 (-03) 


106-115 "• 


14.9 ..75 


0.2-2 it -v- 


W 0.15 


5.11 (-03) 


146-155 


20.3 


0.02 


0.01 


4.65 (-04) 


SS-75 


26.0 


0.12 


0.04 


2.79 (-03) 


206-215 


33.0 


0.01 


0. 004 


2.32 (-04) 



List of experimental data shouing sample averages, standard 
deviations, heights, and velocity conditions measured. 
Density values are calculated by dividing the average sample 
ueight by the sample trap volume. 



Fig. 19 shous the relationship between the average particle 
density caught in the sample trap and the sample trap height above 
the bed surface as a function of Uo/Umf. 



For the relatively low fluidization velocities used during 
the data measurements, (maximum Uo/Umf = 3.81) the freeboard 
height can be assumed infinite. Using this assumption implies 
that all of the particles return to the bed and none are 
elutriated, le complete reflux. Under these conditions, the 
following equation has been suggested by Lewis et al C 3] and Kunii 
and Levenspiel [1] to model the particle loading within the 
freeboard . 
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Particle Qensity (yr»/cn 3 




Freeboard Height (cn) 

Fig. 19 Plot of Particle Density vs Freeboard Height 
as a Function of Uo/Umf . This data was 
Obtained Using the Particle Sampler. 
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Pr 



( 10 ) 



= Po exp{-a 1> 



where: 

Pr = Particle density at height 1 
Po = Particle density at (+)bed surface 
a = Characteristic particle decay length 
1 = Height above bed surface 



Table 5 shows the values for the parameters Po and a obtained 
by linear regression for the curves shown in Fig. 12. 
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TABLE 5 



Uo/Unf 


Po 


a 


Correlation 

Coefficient 


3.81 


0.Z38 


0.1097 


0.379 


3. 17 


0.137 


0.1022 


0.344 


Z . 53 


0.040 


0.1181 


0.337 


Z.3Z 


0.031 


0.1399 


0.857 




Pr = Po 


exp-C-a 1> 




Results of 


linear regression analysis for particle 



loading density (grams/cm ) vs height above the bed 
surface (cn). 



The parameter Po physically represents the particle loading 
density uhich uould be obtained if the sample uere taken at the 
surface of the bed. This is not necessarily the case as is 
indicated by the computer model uhich is discussed in chapter VI, 
but is only a parameter descibmg the particle loading 
distribution in the region of the data obtained. The dashed lines 
belou freeboard heights of 4 cm indicate the region in question. 
The parameter a is a characteristic length of decay for the 
particle flux. 

A correlation betueen the values for Po in table 5 and Uc/L‘mf 
is shoun in Fig. 20. This plot shous that Po is closely related 
with Uo/Umf . Po varies uith (Uo/Umf - 1) approximately to the 2.9 
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Fig. Z0 Plot of Po vs (Uo/Unf -I) Showing Strong 
Dependence of Po on Uo. 
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power, where as Ualsh et al [10] determined the coefficient to be 
about 2.1. Uen and Chen [5] reported a correlation for particle flux 
which is proportional to bubble diameter and (Uo/Umf - 1 > to the 
S/2 power. These correlations are listed in table 6. 

TABLE 6 

- ~ Present work: B 0.13~< Uo/Umf - 1) kg/ m sec 

Uen and Chen B = 1.34 (04) ( Uo-Umf ) kg/ m sec 
Ualsh et al B = 18 ( Uo - Umf ) kg/ m sec 

Comparison of least square fit relations for Po 
as functions of Uo and Umf. 



The differences in these correlations are due to the 
different bed configurations in which the data was taken and the 
measurement technique used. The present work utilized a bed with 
a relatively closely spaced tube configuration and steel grit 
(S.G. 8.1, median size 230 microns) for the bed mass. Both 

ascending and descending particle fluxes were captured in the 
sample. The work of Ualsh et al [10], used a bed with two (2) 
widely spaced horizontal serpentine tubes and Ottawa sand (S.G. 
2.6, median size 755 microns). Also, only descending particle 
flux was used in determining their relations. The correlations of 
Uen and Chen [51 are a result of studies conducted on previous 
research using cylindrical column beds and low mass bed materials 
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(S.G. 0.9 - 2.6). The data for these analysis is based mainly on 

pressure measurements - 

It has been observed by Leuis et al [3] and Uen and Chen £53 
that a is not a strong function of Uo. Both of these studies 
recommend that the characteri stic particle decay length, 1/a, 
could be approximated by an expression of the form: 



1/a = C Uo 



( 11 ) 



Table 7 is a list of correlations obtained by other studies 
and in the present uork. The study by Leuis et al E 31 uas 
conducted uith 75 micron glass spheres in a cylindrical bed. Fig. 
21 shous the relationship between 1/a and Uo in the present uork. 



Present uork 
Leuis et al 
Uen and Chen 
Ualsh et al 

Comparison of 
1/a as functi 



TABLE 7 

1/a = (0.19 
1/a = (1.42 
1/a = (0.25 
1/a = (0.32 

least scuare f 
ns of Uo (m/s) 



0.03 s)Uo m 
0.14 s)Uo m 
0.09 s)Uo m 
0.05 s)Uo m 

t relations for 
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Dependence of 1/a on Uo. 



Particle Size Distribution 



Two methods for determining particle size distribution uere 
used- The sieve method, described in Appendix E, uas used for bed 
material analysis only, due to the small sample sizes obtained from 
the trap. To analyze the small trap samples, an image analyzer 
uas utilized. The procedure used uith the image analyzer is 
described in the experimental procedure chapter. 

The sieve data representing bed mass as a function of 
particle size, is listed in Appendix E. Fig. 22a shows the bed 
mass data converted to particle number as a function of particle 
diameter uhere the interval of particle diameter between 
successive measurements is 10 microns. These plots are faired 
from the data listed in Appendicies E and J. Fig. 22b shous the 
bed mass distribution. The conversion from a mass distribution to 
a particle number distribution uas calculated assuming that each 
particle uas spherical in shape. The volume correspond! ng to a 
given particle diameter uas multiplied by the particle density to 
get a unit particle mass. The mass fraction of the sieve analysis 
corr espondi ng to the specified particle diameter uas then divided 
by the unit particle mass to obtain the represent i t i ve particle 
number. The overall resulting particle number distribution curve 
uas then normalized uith rescect to a maximum value of 100. 

The particle size distribution data obtained from the image 
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Cumulative Mass 



Image Analyzer Analysis 



Sieve Analysis 




~ O 100 200 30O 500 6 00 

Particle Diameter (microns) 

a) Particle Distributions 




b) Mass Distribution 

Fig. 22 Particle Size vs Mass di stri but i on 

and Particle Number of Bed Material. 
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analyzer is listed in Appendix J. The image analyzer uas used to 
analyze the bed mass and sample data from the Uo/Umf = 3.81 data 
set. The bar graphs and the cumulative percentage plots shou the 
number of particles vieued by the image analyzer plotted as a 
function of normalized particle diameter. 

Fig. ZZa shows a plot of the particle number distribution as 
a function of particle diameter as determined by the image 
analyzer. The data is normalized with respect to 100 and is 
compared with the data as determined by the sieve analysis. The 
descrepancy between the two plots can be explained by the methods 
used to determine the respective data. The analyzer first 
determines the cross sectional area of the viewed particle. A 
circle, having the same cross sectional area, is then calculated. 
This results in averaging the smaller minimum diameter with the 

t 

larger maximum diameter of all particles with long cylindrical or 
ellipsoidal shapes. In the sieve however, a large portion of the 
particles will pass through the sieve screen by means of the small 
cross sectional area presented by their longitudinal direction. 

To be consistant in the following sections, the image analyzer 
data will be used to correlate all particle density distributions. 

The particle size distribution as a function of height above 
the bed surface was evaluated using the data collected at Uo/Umf = 
3.81. The image analyser was used to determine the distribution 
by viewing three (3) random samples from each of the six (S) 
height positions. Table 8 lists the various statistical values 
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for the particle size distribution data obtained (Appendix J). 



TABLE 8 



Freeboard 


Average 


fledi an 


Mode 


Diameter at 


Height 


Diameter 


Diameter 


Di ameter 


50Z of 
Cumulative 
f^imber 

Distribution 


( cn) 


( mi crons ) 


( microns ) 


( microns ) 


C microns ) 


3.8 


198 


192 


142 


158 


7.6 


219 


214 


187 


203 


12.3 


204 


189 


162 


167 


17.8 


179 


164 


154 


179 


22.2 


197 


179 


152 


177 


31.1 


212 


197 


194 


190 



Staistical values for particle number distribution 
as a function of freeboard height. A complete 
listing of the data is given in Appendix J. 



Fig. Z3 shous these values plotted against the bed height at 
uhich they uere taken. The excursion of points at the 7.6 cm 
height is assumed to be due to analysis error. The linear 
regression lines for the average and median values show that they 
are ueak functions of collection height. For the particle 
diameter values representing the mode and the 50Z point on the 
cumulative percentage plot, the linear regression lines snow a 
stronger dependence on collection height. 
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Freeboard Height (cm) 

Fig. 23 Analysis of Variation of Particle 
Size vs Freeboard Height. 



Oscilloscope Trace Analysis 



Figs. 24 thru 27 shou examples of the oscilloscope traces 
given in Appendix K. Figs. 24 and 25 represent the traces 
obtained during lou Uo coditions. Under these conditions, the 
normal bubble probe output is high (+5 volts dc ) . When a bubble 
errupts, particles are thrown from the surface of the bed and 
eclipse the light path at the tip of the probe (trace goes to 
zero). Figs. 26 and 27 represent the traces obtained at higher 
Uo conditions when the bubble probe is normally eclipsed (bubble 
probe output is lou, zero) by particles. The presence of a bubble 
is detected by a high output from the bubble probe due to the 
bubble creating a void through uhich the light beam can pass. The 
anemometer output just above the bubble probe is used to determine 
the velocity of the gas jet leaving the erupting bubble. As seen 
in the figures, there is not always a bubble associated with a gas 
jet and vice versa. This is due to the anemometer probe being 
slightly offset from the bubble probe as described in the 
experimental procedure section. 

The gas jet in Fig. 24 is delayed 5 ns and the gas jet in 
Fig. 25 is delayed 35 ms from the point where the bubble traces 
begin their excursion to the zero (eclipsed) condition. Since the 
seperation between the bubble probe and the anemometer is 1 cm, 
this results in an estimated jet velocity of 200 cm/ s (Fig. Z4) 
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Fig. 24 Oscilloscope Trace of Bubble Probe, 
Anemometer Probe, and Solenoid 
=- =~ -=- =’ -^ - Actuation at Lou Uo . 




Fig. ZS Oscilloscope Trace of Bubble Probe, 
Anemometer Probe, and Solenoid 
Actuation at Lou Uo. 
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Fig. Z6 Oscilloscope Trace of Bubble Probe, 
Anenoneter Probe, and Solenoid 
— Actuation at Higher Ho. 




Fig. Z7 Oscilloscope Trace of Bubble Probe, 
Anenonet er Probe, and Solenoid 
Actuation at Higher Uo. 
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and Z8 cn/s (Fig. 25). Using the calibration curve in Appendix 
fl, the velocities are determined to be 518 cm/s for the jet in 
Fig. 24 and 792 cm/ s for the jet in Fig. 25. The large 
difference in these values could be due to a bulge in the bed 
surface from the bubble underneath eclipsing the bubble probe. 

If this is the case, there uould be a delay betueen the time at 
uhich the bubble probe was eclipsed and the eruption of the bubble 
from the bed surface. 

An analysis of Figs. 25 and 27 results in the same 
discrepancies betueen velocities calculated from delay times and 
measured by the anemometer probe. Fig. 27 shows another phenomena 
which occurs quite frequently. The gas jet appears before the 
bubble- is detected. There are several possibilities which can 
explain these occurances. First, the bubble occurs off center 
from the bubble probe. Under this condition, the bubble may erupt 
and initiate a jet which is registered by the anemometer. The 
bubble then continues to rise and the side of the bubble is 
registered by the bubble probe. This explanation can be altered 
to include bubbles coalescing below the surface. Under this 
condition, a bubble may be pulled into the vortex of an already 
eruoting bubble and it eclipses the bubble probe. 

An analysis of the jet velocities and durations for all the 
oscilloscope traces in Appendix K , resulted in the data listed m 
table 9. In later calculations, the values to be used for the jet 
velocity and duration will be 609 . S cm/ s (20 ft/s) and 20 ms 
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respectively. 



TABLE 9 







Average 


Stnd 


Dev 


Jet 


Velocity 


689 


cm/s 


533 


cm/s 






22.5 


ft/s 


18. 


ft/s 


Jet 


duration 


21.25 


ms 


13.22 


ms 



Average and standard deviation of jet velocity 
determined from oscilloscope traces in 
Appendix K. 



Sample Ueioht Versus Bed Activity Correlation 

The main purpose of inserting the bubble and anemometer probe 
beneath the sampling apparatus- was to determine whether a 
correlation exists between the sample weight collected and the 
presence of bubble eruptions and gas jets. To evaluate the 
photographs taken of the oscilloscope traces (Appendix K), the 
following information was required- 

A. The average velocity of the particles as they 
travel from the bed surface to the trap. 

B. The distance the particle must travel to reach 
the trap. 

C. The closure time of the trap relative to the 
bubble eruptions and -gas jets. 



89 



h e average particle velocity uas calculated by deternming 
tne average particle size and using the height output from the 
comout er trajectory model. The average particle size as 
aetemined by the image analyzer uas 180 rnicons. To determine the 
average velocity using the model, the following initial conditions 
were input to the program: 

Superficial velocity ( Uo ) = 57.9 cm/s (1.9 ft/s) 

This uas the actual velocity 
measured during sampling. 

Jet velocity ( Llj ) = 609. S cm/s (Z0 ft/s) 

This value uas determined from 
the anemometer data. 

Jet duration (t jet) = 20 ms 

This value uas determined from 
the oscilloscope traces. 

The height attained by a 180 micron particle, as determined 
by the model uith the above conditions is 44.7 cm in 0.295 
seconds. This results in an average velocity of 44.7/0.295 - 
151.5 cm/ s (4.97 ft/s). 

The distance a particle must travel to reach the center of 
the t rap from the bed surface is obtained from the data listed in 
Appendix I. Since trao height is measured from the bed surface to 
the bottom of the trap, 4 cm (0.13 ft) must be added to the trap 
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heights to obtain the distance to the center of the trap. The 
».otal tine required for the particles to leave the bed surface and 
arrive at the center of the trap can nou be deternined. Table 10 
shous the transit tine for each set of data for which the 
oscilloscope traces were photographed (Appendix K). 



TABLE 10 



Data 


Trap 


Trap 


Total 


Tine 


set 


Height 


Height 


Transit 


Bef ore 


Mjmber 


3otton 


Center 


Tine 


Actuation 




(cm) 


( cm ) 


( ns ) 


( ns ) 


56-65 


24.15 


27.13 


179 


136 


66-75 


26.03 


30 . 03 


198 


155 


76-85 


22.23 


26.23 


173 


130 


86-95 


24.76 


28.76 


190 


147 


96-105 


13.65 


17.65 


116 


73 


List of 


transit tines for 


particles traveling from 


the bed 


surface to 


the center of the trap 


. The tine 


prior to actuation 


of the 


sanple trap is 


also shown 


(Total 


tine - 42.6 


ns ) . 







The total closure tine of the sampling apparatus (closure of 
solenoid power supply switch SI until trap is shut) is deternined 
in Appendix D to be 42.6 ns. This value is subtracted from the 
particle transit tine to obtain the tine before actuation in which 
particles leaving the bed surface will be caught in the trap. 

To deternine whether or not a correlation exists between the 
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weight of a sample and the amount of bed activity present in the 
bed, a weighted analysis was used. Each sample was evaluated 
three ways. First, if the sample weight was less than the average 
sample weight, a weight factor (U) of -1 was assigned to it. If 
the sample weight was greater than the average, a weight factor of 
+1 was assigned. Second, the oscilloscope trace was analyzed at 
the point corresponding to the "time before actuation*' listed in 
Table 10. This position is the tine before the solenoid power 
supply trace is present (Fig. 28). If the trace showed signs of 
a bubble eruption or a gas jet at this point, the position factor 
(p) was assigned a value of +1. If no activity was present, the 
position factor was assigned a value of -1. Third, each trace was 
analyzed again at the "before actuation tine" but bed activity 
within a +/- 25 ns region was counted. If there was bed activity 
within this region, the region factor (r) was assigned the value 
of +1. If no activity was present, the region factor was assigned 
the value of -1. The following equation was used to assign a 
correlation factor to each sample- 

Q=Up+Ur (10) 

where - 

Q - Correlation factor 

U = +1 if sample is > average of data set 

-1 if sample is ( average of data set 

p = -1 if no activity is present in bed 

+1 if bubble or gas jet activity present 



Tine ( ns ) 

0 100 Z00 300 400 500 




Fig. Z 8 Typical Oscilloscope Trace Qjnng 
Sanpling Procedure. 
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r = “ 1 if no activity is present +/- 25 ns 
+1 if activity is present +/- 25 ns 



Table 11 shows the correlation factor for each sample. The 
resulting average correlation factor for all samples is 0.32, 

This value suggests that no correlation can be made between the 
samole weight and bed activity. There are several explanations 
for this. Bubbles not detected by the bubble probe or debris from 
previous bubbles returning to the bed surface can influence the 
sanple weight by increasing the amount of particles collected, 
Also, bubbles of smaller size nave a smaller velocity and would 
therefore produce a particle stream which is either slower or just 
delayed in leaving the bed surface. This would result in lower 
sample weights than expected. 



94 



TABLE ii 



Sampl e 
Number 

So 

57 

S3 

S3 

50 

SZ 

53 

54 

So 

S7 

S3 

70 

77 

73 

83 

84 

8G 

88 

83 

30 

9Z 

93 

94 

95 

96 

97 

98 

99 
1C0 
101 
10Z 



U 

1 

-1 

1 

1 

-i 

1 

-1 

-1 

1 

-1 

-1 

-I 

1 

1 

-1 

1 

-1 

-1 

-1 

1 

-1 

1 

•1 

-1 

1 

1 

1 

1 

-1 



P 

1 

1 

1 

- 1 
-i 
1 
1 

-1 

-1 

1 

1 

1 

1 

1 

1 

1 

-1 

1 

1 

-1 

-1 

1 

-1 

1 

-1 

-1 

-1 

-1 

-1 

1 

-1 



r 

1 

1 

1 

-1 

-1 

1 

1 

-1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



-1 

1 

-1 

-1 

1 

1 



Q 

Z 

-Z 

z 

Z 

-z 

-7 

z 

z 

0 

z 

-z 

-z 

z 

-z 

z 

z 

0 

z 

-z 



z 

0 

z 

0 

z 

0 

-z 

-z 

z 

0 



Average = 0. 3Z 



List of samples and their correlation parameters. The 
resulting average value for Q indicates that no 
correlation can be made from the data obtained to 
indicate by sample ueight whether or not any bed 
activity occured belou the sample trap. 
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CHAPTER VI 



TRAJECTORY flQDFI RESULTS AND DISCUSSION 
Selection of Baseline Parameters 

The input to the model consists of the following five (5) 
parameters - 



1) Superficial bed velocity ( Uo ) 

Z) Initial particle velocity ( Upo ) 

3) Peak gas jet velocity ( Uj ) 

4) Gas jet duration (tj) 

5) Particle distribution of the bed mass 

The baseline values for each of these inputs was determined 
to be as close to the actual experimental bed conditions as 
possible. 

Superficial Bed Velocity 

The superficial bed velocity uas determined directly from the 
experimental data. For the bed conditions discussed in this 
section, the superficial velocity used is the same as the highest 
velocity condition under uhich the particle sampler uas used. 

The superficial velocity within the bed for these samples was 
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calculated to be 57.3 cn/s (1.3 ft/s). 



Initial Particle Velocity 

To determine the initial velocity of the particle, it uas 
assumed that the particle uas located at the nose of a bubble and 
would therefore have the bubbles velocity. To determine the 
bubble velocity, Kunii and Levenspiel [1] give the following 
equat ion- 



Ub = 0.711 Cg Db] + Uo - Umf 



( 11 ) 



□here- 



Ub = Bubble velocity 
Uo - Superficial bed velocity 
Umf = Minimum fluidization velocity 
Db = Bubble diameter 
-g = gravi atat ional acceleration 



Observations of the bed material cLiring the sampling 
operation suggests that the average bubble diameter present in the 
bed is approximately S cm (Z.4 in). For Uo = 57.3 cm/s, Umf = 

15.2 cm/s and Db = Sen, the bubble velocity is calculated to be 

37.2 cm/s <3.13 ft/s). This value uas used as the initial 
particle velocity for the base line data. 
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Peak Gas Jet Velocity 



To determine an average peak jet velocity, the anemometer 
output on the oscilloscope traces (Appendix K) were analyzed- To 
determine the relationship between the voltage output from the 
anemometer (which is displayed on the oscilloscope) and the actual 
gas jet velocity, the anemometer was calibrated in a wind tunnel - 
The resulting calibration curve for the anemometer output, as 
displayed on the oscilloscope trace, is given in Appendix tt. The 
average peak velocity of the anemometer traces analyzed in chapter 
V was determined to be approximately 609.S cm/s (20 ft/s). 

Gas Jet Duration 

To determine an average gasrjet duration time to use as a 
baseline, the anemometer traces in Appendix K were used. The time 
duration of each gas jet analyzed is determined directly from the 
oscilloscope trace. The average duration of a gas jet was 
determined to be approximately 20 ms in chapter V. 

Particle Distribution of the Bed Mass 

Two particle distributions for the bed mass were available 
for use as the baseline data. Both particle distributions were 
determined using the same sample but analyzed using different 
analysis techniques- The two methods used are the sieve method 
and the image analyzer method, both of which are described earlier 
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in chapter V. The tuo distributions are shown in Fig. 29. 

The bed distridu tion determined by the image analyzer was used to 
be consistent with the particle distributions measured at various 
freeboard heights. 



Typical Output Using Baseline Parameters 

Fig. 30 shows the resulting maxinun particle height obtained 
as a function of particle diameter for the baseline conditions. 

The general shape of the curve is deternined by the dominating 
force acting on the particle. For large particles, the dominating 
force at these air velocities is gravity. Therefore, momentum 
(initial particle velocity) is the controlling factor determining 
the maximum height attained by the particle. As the particle size 
decreases, the proportion of drag force to gravitational force 
becomes larger. For the smaller particles, the drag force becomes 
the dominate force. The result is that the maximum height 
attained by a particle continues to increase as particle size 
decreases. This continues until the value of Uo approaches the 
terminal velocity of the smaller particles. At this point, the 
particle is totally dominated by the air flow in the bed. This is 
seen in the small decrease in height attained by particles less 
than 110 microns as a result of the particles rapid deaccel erat ion 
tc Uo after the gas jet has stooped while the larger particles 
continue a little higher due to momentum. The particle height 
continues to decrease for a short time until Uo becomes equal to 
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Rul alive Pcirliclti Number 



red Data Front Appendix E and J. Interval of 10 microns. 




Fig. 29 Relative Particle Number Distributions 
of Bed Material by Sieve and Image 
Analyzer Analysis. 




Particle Qia. (Microns) 



Fig. 30 



Maximum Particle Height vs Particle 
Diameter for Baseline Conditions. 



or greater than the terminal velocity for the remaining particles, 
hese particles would be elutriated from the system. 

Figs. 31 thru 36 show the individual particle densities at 
different heights above the bed. All of the heights shou an 
increase in particle density for the smaller diameter particles. 
This is the result of the smaller particles falling at their 
respecTive terminal velocities uhich is slower than for the larger 
particles. the result is that the smaller the particle velocity, 
the longer the particle exists within a given height region and 
thus, the larger the particle density. It is also due to the 
larger number of smaller particles present within the system. 

Fig. 37 shows the particle density distribution above the 
bed. The two curves represent the effect of varying the range of 
particle- diameters used in the model- The two ranges are 80 * 570 
microns and 80 - 1070 microns. The peak of the 80 - 570 micron 
distribution is at about 18 cm of bed height whereas the other 
distribution peaks at about 11 cm. The region shown on the graph 
below these heights is the area known as the splash zone. The 
shape of this part of the curve is due to the initial acceleration 
of the particles leaving the surface of the bed by the gas jet. 

The particles then begin to slow down due to drag and gravity 
resulting in an increase in particle density. The peak on this 
curve coincides with the maximum height attained by the largest 
particles, and therefore, it can be assumed that the change in 
slope is due to the loss of particles as they return to the bed. 



101 



100 



fe 30 - 

n 

5 30 - 

2 

a 73 -■ 




Fig. 31 Relative Particle Number vs Particle 
Diameter for Baseline Conditions. 

Freeboard Height of 4 cm. 
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Fig. 3Z Relative Particle Nunber vs Particle 
Oi anet er for Baseline Conditions. 
Freeboard Height of 8 cn. 
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Fig. 33 Relative Particle Number vs Particle 
Diameter for Baseline Conditions. 
Freeboard Height of 12 cm. 
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Fig. 34 Relative Particle Number vs Particle 
Diameter for Baseline Conditions. 
Freeboard Height of 18 cm. 
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Fig. 35 Relative Particle Number vs Particle 
Diameter for Baseline Conditions. 
Freeboard Height of 22 cm. 
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Fig. 36 



Relative Particle Numoer vs Particle 
Diameter for Baseline Conditions. 
Freeboard Height of 31 cm. 
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Conditions . 




Fig. 38 In Particle Density/Unit Volume 
vs Freenoard height for Baseline 
Condi t ions . 
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Fig. 38 is a semi-ln plot of the sane data as plotted in Fig. 

37. The sharp drop at the right end of the plot is due to 
analyzing particles with a nininun diameter of 80 microns. Had 
smaller particles been included, the slope would have approached 
zero instead of infinity which would model the elutriation of 
particles. The slope of the line to the right of the peaks shown 
in Fig. 38, which are located at 11 and 18 cm of height above the 
bed surface are the same and are approximately -0.117 grams/cm . 

To decrease computation time, the 80 - 570 micron particle 
distribution was used in the following sensitivity analysis as the 
decay slopes are equal. 

Model Sensitivity Analysis 

An analysis of the change in bed characteristics due to a 
change in one of the f ive~model input' parameters was conducted. 

The parameters changed and the values used are as follows* 



Superficial velocity <Uo) : 

Baseline Uo : 57.9 cm/ s (1.9 ft/s) 



1. 30.5 cm/s ( 1 ft/s) 

2. 91.4 cm/s (3 ft/s) 



Initial particle velocity ( Upo ) : 

Baseline Upo- 



1 . 
2 . 



119.2 cm/s (3.19 ft/s) 

61 . 0 cm/s (2 f t/s ) 
152.4 cm/s ( 5 ft/s ) 



1 06 



3. 305.0 cn/s (10 ft/s) 



Gas Jet velocity ( U j ) : 

Baseline Uj - G09.S cn/s (20 ft/s) 



1. 305.0 cn/s (10 ft/s) 

2. 457. Z cn/s (15 ft/s) 



Gas jet duration (tj): 

Baseline t j : 20 ns 



1. 10 ns 

2. 30 ns 



Particle distribution of bed mass- 

Baseline distribution- Inage Analyzer 



1. Sieve analysis data 



The bed character ist i cs uhich are evaluated- against the 
baseline character istics are: 



1. The naxinun height attained by each particle from 80 to 
570 nicrons. 

Z. The density of each particle size (80 - 570 nicrons) at 
freeboard heights of: 



a. 4 cn (1.5 in) 

b. 8 cn (3.1 in) 
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c. 12 cm (4.7 in) 

d. 18 cm (7.1 in) 

e. ZZ cm (8.7 in) 

f. 31 cm (1Z.Z in) 

3. The particle density distribution in the freeboard. 



Variation of Superficial Velocity (Uo) 

Fig. 39 shous the influence of Uo on the maximum height 
attained by the particles. As uas described in the previous 
section, the larger particles are dominated by momentum and not 
drag. ZJhis is shown in Fig. 32 by the small change in maximum 
height attained by the large particles due to a change in Uo. As 
the particle size decreases, "the drag force becomes the dominate 
force and the effect of Uo on particle height increases. As seen 
in Fig. 39, a change in Uo produces a moderate change in the 
maximum height attained by the small particles. 

Figs. 40 thru 45 show the effect of Uo on individual 
particle densities at different heights above the bed. All of the 
heights show an increase in particle density for the smaller 
diameter particles as Uo is increased- This is the result of the 
smaller particles falling at their respective relative terminal 
velocities which is slower for smaller particles. A decrease in 
the maximum particle diameter present at a given freeboard height 
as Uo is decreased can be observed in Figs. 40 thru 45- This 
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FiQ. 39 Maxinun Particle Height vs Particle 
Diameter as a Function of Uo. 




Fig. 40 Relative Particle Number vs Particle 
Diameter as a Function of Uo. 
Freeboard Height of 4 cm. 
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Diameter as a Function of Uo. 
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Diameter as a Function of Uo. 
Freeboard Height of 12 cm. 
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Relative Particle Number _ Relative Particle Number 



100 -■ 
30 - 



91.4 cn/s 



30 - 
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57.9 cn/s 
30.5 cn/s 




Fig. 43 Relative Particle Nunber vs Particle 
Dianeter as a Function of Uo. 
-Freeboard Height of 18 cn. 




Fig. 44 Relative Particle Nunber vs Particle 
Dianeter as a Function of Uo. 
Freeboard Height of 22 cn. 
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Fig. 45 Relative Particle Number vs Particle 
Diameter as a Function of Uo. 
Freeboard Height of 31 cm. 
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decrease is due to the decrease in lift given to the ascending 
particle as Uo is decreased. As a result, all particles achieve a 
louer maximum height uhen Uo is decreased. 

Fig. 4G shous the effect of Uo on the particle density 
distribution above the bed. There is a slight increase in overall 
particle density uith an increase in Uo. As expected, there is 
also an increase in the naxirum height attained by the particles 
uith an increase in Uo . Fig. 47 is a semi-ln plot of the same 
data plotted in Fig. 46. Table 12 lists the slopes of the lines 
to the right of the peak shoun in Fig. 47 unich is located at 
about 18 cm of height above the bed surface. This table shous an 

y 

increase in the particle distribution slope as Uo is decreased. 



TABLE 12 



Uo 

cm/s 



Slope 

gms/cm 



I ntercept 
gms/ cm 



30. S 
57.9 
31.4 



-0.136 

-0.117 

-0.098 



7.036 

6.758 

6.473 



Effect of Uo on the slope of the 
particle density distribution as a 
function of height for the 
distributions shoun in Fig. 47. 
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In Particle Oensi ty/lini t Vul 




Fig. 46 Particle Density/Unit Volume vs 
Freeboard height as a Function 
of Uo. 




Fj>Q. ** s In Par* i ci e Density /Unit Volume 

vs Freeboard height as a Function 
of Uo. 
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Variation of Initial Particle Velocity ( Upo ) 



Fig. 48 shows the relation between initial particle velocity 
(Upo) and the maxinun height attained by particles of different 
size. This figure also shows that the larger particles are mainly 
momentum dependent while the smaller particles are drag dependent. 
If Upo were increased further, the larger particles would continue 
to increase their maximum height. The smaller particles would 
approach a height limit which is dependent upon the superficial 
velocity in the bed. 

Figs. 49 thru 54 show the effects of increasing Upo on the 
individual particle densities for increasing heights. In general, 
the individual particle distributions undergo the same relative 
changes from the bed surface to the maximum height position. The 
difference being the height above the bed surface at which the 
particular distribution is present. 

Fig. 55 shows the effect of varying Upo on the particle 
density distribution in the freeboard. As was noted previously, 
the peak density occurs at the point where the largest particles 
attain their maximum height above the bed. As Upo is increased, 
the curve to the right of the peak becomes shorter and steeper. 
This trend will continue until the the effect of the smaller 
particles returning to the bed cause the slope to decrease. Fig. 
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Freeboard Height (cn) 




Fig. 48 Maxinun Particle Height vs Particle 
Oianeter as a Function of Upo. 
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Fig. 49 Relative Particle Number vs Particle 
Diameter as a Function of Upo. 
Freeboard Height of 4 cm. 




Fig. 50 



Relative Particle Number vs Particle 
Diameter as a Function of Upo. 
Freeboard Height of 8 cm. 
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Relative Particle Number , Relative Particle Number 
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Fig. 51 



Relative Particle Number vs Particle 
Diameter as a Function of Upo. 
Freeboard Height of 12 cm. 




Fig- 52 Relative Particle Number vs Particle 
Diameter as a Function of Upo. 
Freeboard Height of 18 cm. 
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Fig. 53 Relative Particle Number vs Particle 
Diameter as a Function of Upo. 

___dr=^ Freeboard Height of 22 cm. 
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Relative Particle Number vs Particle 
Diameter as a Function of Upo. 
Freeboard Height of 31 cm. 
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5o is a seni'ln plot of the data in Fig. 55. Table 13 confirms 
what the slope is becoming steeper as Upo is increased although 
the influence of initial particle velocity on the slope of the 
density vs height curve is moderate. 



TABLE 13 



Upo 


Slope 


Intercept 


cm/s 


gms/cn 


gms/cm 


61.0 


-0.118 


6.058 


97. Z 


-0.117 


6.758 


15Z.4 


-0.1Z5 


8.500 


304.8 


-0.Z00 


18. Z0 



Effect of Upo on the slope of the 
— particle density distribution as a 

function of height for the 
distributions shown in Fig. 56. 



Variation of Jet Velocity ( U j ) 

Fig. 57 shous the effect on maximum particle height uhen Uj 
is varied. As in the tuo previous analysis’, the influence of 
momentum and drag on the different particle sizes is apparent. A 
decrease in jet velocity from the baseline value results in a 
large drop in small particle height and a small change in large 
particle height. This suggests that jet velocity has a strong 
effect on the ciscersion, or seperation of particles of different 
sizes at increasing freeboard heights. Increasing Uj results in a 
stretching effect of the particle size distributions above the 
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